o 

X 



Transient Growth in Shear Flows: 
Linearity vs Nonlinearity 

Chris C.T. PringleQ and Rich R. Kerswelfj] 

Department of Mathematics, University of Bristol, 

University Walk, Bristol BS8 1TW, United Kingdom 

(Dated: May 11, 2010) 

Two approaches to the problem of transition to turbulence of shear flows are popular in the lit- 
erature. The first is the linear one of transient growth which focuses on the likely form of the most 
'dangerous' (lowest energy) turbulence-triggering disturbances. The second is the nonlinear calcula- 
tion of the laminar-turbulent boundary which instead focuses on their typical amplitudes. We look 
to bridge the gap between these two perspectives by considering the fully nonlinear transient growth 
problem to estimate both the form and amplitude of the most dangerous disturbance. We thereby 
discover a new nonlinear optimal disturbance which outgrows the well-known linear optimal for the 
same initial energy and is crucially much more efficient in triggering turbulence. The conclusion is 
then that the most dangerous disturbance can differ markedly from what traditional linear transient 
growth analysis predicts. 
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Shear flows are ubiquitous in nature and engineering, 
and understanding how and why they become turbulent 
has huge economic implications. This has led to a num- 
ber of simplified canonical problems being studied such 
as plane Couette flow, channel flow and pipe flow which 
commonly exhibit turbulent behavior even when the un- 
derlying laminar state is linearly stable. In this case, a 
finite amplitude perturbation is required in order to trig- 
ger turbulence and a leading question is then what is the 
'most dangerous' or 'smallest' such perturbation (with 
the metric typically being energy). Beyond its intrinsic 
interest, such information is fundamentally important for 
devising effective control strategies to delay the onset of 
turbulence. 

Historically, linear transient growth analysis [lMa] has 
been used to identify dangerous disturbances which are 
efficient at triggering turbulence. This focuses on a linear 
mechanism whereby infinitesimally small perturbations 
can interact with the underlying shear profile in order 
to create much more energetic disturbances. Ultimately, 
these disturbances ebb away if the shear profile is linearly 
stable but the thinking is that the large growth possible 
can catapult the disturbance into a regime where non- 
linear effects sustain its energy away from zero for all 
times. 

Alternatively, recent progress has been made [0, Lfl] 
in numerically tracking the laminar-turbulent boundary 
which represents the surface in phase space which sepa- 
rates those initial conditions which will trigger a turbu- 
lent episode from those which will simply relaminarise. 
Since the tracking technique hinges upon carefully se- 
lecting initial conditions and integrating forward in time, 
the part of the boundary revealed is effectively confined 
to the neighborhood of the limiting set of the boundary- 
confined flow dynamics. This is found to have signifi- 
cantly higher energy levels than that actually needed to 
trigger turbulence by carefully-tuned initial disturbances 



(e.g. [8J). This disparity in energy levels is none other 
than an expression of the large transient growth endemic 
in shear flows but now crucially translated into the non- 
linear regime. The most dangerous disturbance corre- 
sponds with the minimum energy point on the laminar- 
turbulent boundary and stands to gain the largest energy 
as it sweeps up to the limiting set energy plateau. In this 
Letter, we pursue the promising strategy of extending 
the usual linear transient growth analysis into the non- 
linear regime to identify the most dangerous disturbance 
in pipe flow. Our findings should be equally relevant to 
other shearing flows such as plane Couette flow, channel 
flow, boundary layers, etc. 

The transient growth problem is the optimisation ques- 
tion: what initial condition u(x, t = 0) (added as a per- 
turbation to the laminar flow 2(1 — 4s 2 )z) for the gov- 
erning Navier-Stokes equations with fixed (perturbation) 
kinetic energy Eq will give rise to the largest subsequent 
energy Et at a time t = T later. This corresponds to 
maximising the functional 

X := <±u(x,T) 2 )-A<±u(x,0) 2 -£ ) 
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where ( } represents volume integration; (s, <p, z) are 
cylindrical coordinates directed along the pipe; A, i/(x, t), 
II(x, t) and r(i) are Lagrange multipliers imposing the 
constraints of initial energy Eq, that the Navier-Stokes 
equations hold over t E [0,T], incomprcssibility and con- 
stant mass flux in time respectively (the system has been 
non-dimensionalised by the pipe diameter D and the bulk 



velocity U with Re := pUD/p where p is the density 
and p is the dynamic viscosity). Vanishing of the varia- 
tional derivatives requires that u must evolve according 
to the Navier-Stokes equations, v evolves according to 
the adjoint-Navier- Stokes equations and at times t = 
and T we have optimality and compatibility conditions 
linking the two sets of variables (e.g. see [9J for details of 
the linearised problem). The method of solution is one 
of iteration as follows: 

• Make an initial guess for u(x, t = 0). 

• Allow u(x, t) to evolve according to the Navier- 
Stokes equations until t = T . 

• Solve the compatibility condition for i/(x,T), 
<LS?/£u(x, T) = u(x, T) - i/(x, T) = 

• Allow the incompressible field f(x, t) to evolve 
backwards in time until t = via the adjoint- 
Navier-Stokes equations 

dv . o, du 1, . - _. 
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FIG. 1: The evolution of the linear and nonlinear optimals at 
Re — 1750. The blue line corresponds to the nonlinear opti- 
mal for Eq = 2 x 10 -5 while the red line is the linear optimal 
(Eq — > 0). The nonlinear result produces more growth and 
actually reaches its maximum at a slightly earlier time than 
T. 
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• Move u(x, 0) in the direction of the variational 
derivative <Lzf/<5u(x, 0) = — Au(x, 0) +i^(x, 0) and 
repeat. 

Both direct and adjoint equations were solved using 
a fully spectral, primitive variables approach. Time 
stepping was done using a second order fractional step 
scheme, checked carefully against the code of [lfj. The 
computational domain was a short periodic domain of 
length 7T radii with typical spatial resolution of 29 real 
Fourier modes azimuthally, 11 real Fourier modes axially 
and 25 modified Chebyshev polynomials radially in each 
of the 8 physical scalar fields (u, v, w,p, V\, 1/2, V3, II). All 
results have been checked for robustness to resolution 
changes. Retention of the nonlinear terms poses a fresh 
technical challenge: although the adjoint equation is lin- 
ear in v, it is dependent on the evolution history of the 
forward variable u which now must be stored. 

The linear transient growth optimal U; ira (x;_Re) in 
pipe flow is well-known to be streamwise-independent 
(2D) rolls which evolve into much larger streamwise- 
independent streaks [5J: see figs 1 and 2. Maximum 
growth occurs at T Hn w 12.2 x Re/1000 (D/u) [lj. In- 
troducing nonlinearity (ie increasing Eq from 0), setting 
T = 12.2 x i?e/1000 (D/U) and allowing only 2D flows, 
leads smoothly to a modified 2D optimal U2D (x; [Eg , Re) 
with monotonically decreasing growth (see [l2j for an 
equivalent result in boundary layers). By Eq ~ 10~ 3 (in 
units of the laminar flow's kinetic energy, used hence- 
forth), the energy magnification falls to about 50% of 
its linear value at Re = 1750, and continues to decrease 
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FIG. 2: Three snapshots of the linear optimal (top) and six 
snapshots (middle & bottom) of the 3D optimal for Re — 1750 
and Eq — 2 x 10 _J during its evolution. Labels refer to figure 
[T] arrows indicate cross-sectional velocities and colours axial 
velocity beyond the laminar flow (white/light for positive and 
red/dark for negative: outside shade represents zero). 
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FIG. 3: Re = 1750, E = 2 x 10~ 5 . The iterations are seeded 
with a noisy version of the 2D optimal and, after approaching 
two 'saddle' points, eventually converge onto the 3D optimal. 



thereafter. Opening the optimisation up to fully 3D flows 
initially just recovers the 2D result but once E crosses 
a small threshold E 3 £> (w 10~ 5 at Re = 1750), a com- 
pletely new optimal U3£>(x; Eq, Re) appears. This new 
3D optimal emerges from the optimisation procedure af- 
ter it initially appears to converge to the 2D optimal and 
then transiently visits an intermediate state: see fig 3. 
Identifying this 'loss of stability' of the 2D optimal pro- 
vided an efficient way to compute Ezr>{Re). All optimisa- 
tion results were robust over three very different choices 
of starting flow: a) uu n with noise; b) the asymmetric 
travelling wave |13| : and c) a turbulent flow snapshot 
(all rescaled to the appropriate initial energy). 

Given the intensity of the runs (O(200) iterations and 
each iteration requires integrating forwards and back- 
wards over the period [0, T]), we concentrated on two 
values of the Reynolds number, Re = 1750 and 2250, and 
the corresponding energy ranges \E%d = 1-35 x 10~ 5 , 2 x 
10~ 5 ] and [E 3D = 4.8 x 10~ 6 ,6.25 x 10~ 6 ]. Over both 
intervals the 3D optimal has essentially the same appear- 
ance: see fig. HJ Unlike the linear optimal which is glob- 
ally simple in form and undergoes an evolution that is 
well established (rolls advecting the mean shear to gen- 
erate streaks), the 3D optimal is localised to one side of 
the pipe and initially has both rolls and streaks of com- 
parable amplitude. Figures Q] and [2] show a new 2-stage 
evolution: a preliminary phase when the flow delocalises 
followed by a longer growth phase where the flow struc- 
ture stabilises to essentially two large-scale slow streaks 
sandwiching one fast streak near the boundary. 

If Eq is increased beyond the ranges quoted above, the 
iterative procedure fails to converge for reasons which 
are unclear. One possibility is that at these energy lev- 
els, the laminar-turbulent boundary has been crossed and 




FIG. 4: Re=2500. The green line shows the evolution of the 
3D optimal when given initial energy S c . Because it is on the 
laminar-turbulent boundary two nearly identical initial condi- 
tions diverge after, in this case, 220D/U. The blue line is the 
evolution of the 2D optimal for the exact initial energy Sl m 
for which the streaks become linearly unstable. The red line 
shows the 2D optimal given initial energy <? c hn and allowed 
to evolve until it reaches a maximum amplitude whereupon 
0.1% by amplitude unstable perturbation is added. Again the 
laminar-turbulent boundary can be identified. 




FIG. 5: The effect of the initial energy on the growth of Auu n 
(red) and Au 3D (blue) at Re - 1750. For small E the 2D 
result is the optimal but after E^d — 1-35 x 10 , the 3D 
optimal takes over. The vertical dashed line corresponds to 
S c , with the dotted lines being the relevant errorbars. Inset: 
The dependence of Esd (red), S c (blue) and <£ l s ln (green) on 
Re {S^ n is even higher). The error bars shown on S c at Re = 
1750 and 1875 are because the laminar-turbulent boundary is 
difficult to identify at these Re due to the pipe shortness. 



the lack of convergence is due to the flow becoming tur- 
bulent. The ensuing sensitivity to noise would make 
the optimisation non-smooth. However, direct numer- 
ical simulation starting with the 3D optimal does not 
reveal a turbulent episode implying that there is still 
an energy gap between where the 3D optimal emerges 
as the solution to the transient growth problem and the 
lowest ener gy o f any initial condition which can trigger 



turbulence la. What can be shown, however, is that 



the 3D optimal is much more efficient at triggering tur- 
bulence than the linear optimal when rescaled. Taking 
the initial condition Au 3 d(x;2 x 10~ 5 , 1750), we gradu- 
ally increase the rescaling factor A until a critical energy 
Eq = S c {Re) is reached at which turbulence is triggered. 
Calculating the corresponding quantity for the linear op- 
timal turns out to be less clearly defined because some 
3D noise is needed to trigger turbulence. As a result we 
make 2 different estimates, one strictly conservative and 
the other more realistic. The first S' l s m is obtained by 
taking Auu n (x; 1750) and finding the initial energy for 
which the resultant streaks are just linearly unstable in 
this periodic domain [14J. In the second S^ n , the same 
initial condition was used but 0.1% of the most unstable 
perturbation (as found from the previous computation) 
is added to the streaks when they reach maximum am- 
plitude. S\ m should be a (low) conservative estimate but 
even this is O(10) times larger than S c at Re — 2500 - 
see figure [4]- whereas the more realistic # c is O(100) 
times larger. 

In figure [Sunset) we plot Ej,d, $c and <§\ ln as a 
function of Re which emphasizes that the 2D optimal 
(for which the linear result is an excellent approxima- 
tion) ceases to be a global maximum at an energy (at 
least) several orders of magnitude before it approaches 
the laminar-turbulent boundary. The 3D optimal, in 
contrast, crosses the laminar-turbulent boundary only 
shortly after it emerges at E%£>, This indicates that the 
3D optimal provides a rapid means of bridging the gap 
between when the 3D nonlinear optimal surpasses the 
linear result and when turbulence can be triggered. To 
achieve this, the energy growth experienced by the 3D op- 
timal must increase dramatically with Eq which is illus- 
trated in figure [5j It is worth remarking that the lowest 
possible energy to trigger turbulence must be bounded 
below by E 3 d [15] and above by S c . 

In this letter we have demonstrated that including non- 
linearities in the problem of transient growth critically 
changes the result close to the onset of turbulence. Al- 
though we have not been able to calculate these solu- 
tions all the way up to the laminar-turbulent boundary 
as originally hoped, we provide evidence that they are 
very efficient at triggering turbulent episodes, notably 
more so than the linear result. Admittedly, we have only 
considered a short periodic domain and so the natural 
question is what will happen in larger domains. Here 
we expect further localisation of the optimal since en- 



ergy is defined as a global quantity whereas nonlinearity 
is important whereever the velocity field is locally large. 
This strongly suggests that in a long pipe the optimal 
should localise in the axial direction as well (the 3D op- 
timal found here is already localised in the radial and 
azimuthal directions). This squares well with the exper- 
imental observation that small local perturbations can 
trigger high energy global turbulence. More importantly, 
the fact that localised flow structures should emerge from 
this type of theoretical analysis bodes well for a greater 
connection between theory and experiments which natu- 
rally introduce localised disturbances. 

Finding fully nonlinear optimals is a time-consuming 
pursuit due to the slow convergence of the iterative pro- 
cedure and the need to look within small energy windows. 
Their discovery has had to wait almost two decades after 
the linear result was established in pipe flow. As com- 
puter power steadily increases, we envision that these 
new nonlinear optimals will start to come within easy 
reach. That the optimal found here represents some- 
thing entirely different to the previously known linear re- 
sults suggests that they will open up a whole new means 
of triggering turbulence and a whole new way of under- 
standing how transition occurs. 

The calculations in this paper were carried out at 
the Advanced Computing Research Centre, University 
of Bristol. 
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